Cooper pairing and single particle properties of trapped Fermi gases 
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We calculate the elementary excitations and pairing of 
a trapped atomic Fermi gas in the superfluid phase. The 
level spectra and pairing gaps undergo several transitions as 
the strength of the interactions between and the number of 
atoms are varied. For weak interactions, the Cooper pairs 
are formed between particles residing in the same harmonic 
oscillator shell. In this regime, the nature of the paired state 
is shown to depend critically on the position of the chemical 
potential relative to the harmonic oscillator shells and on the 
size of the mean field. For stronger interactions, we find a 
region where pairing occur between time-reversed harmonic 
oscillator states in different shells also. 



I. INTRODUCTION 

The trapping and cooling of trapped fermionic atomic 
gases is becoming a very active area of experimental re- 
search. Presently, temperatures as low as ~ Q.2Tp have 
been obtained for ^^K |] and ^Li §-|l with Tp denoting 
the Fermi temperature. Part of the motivation for this 
impressive experimental progress is to observe the theo- 
retically predicted transition to a superfluid state below a 
critical temperature Tc- With a large and attractive low 
energy effective interaction between two hyperfine states 
trapped, Tc for such a transition should be experimen- 
tally obtainable ||]. 

The purpose of the present paper is examine the el- 
ementary excitations of the ground state for an atomic 
superfluid gas in a spherical harmonic trap. Varying the 
interaction strength and the number of particles in the 
trap a number of interesting transitions in the level spec- 
tra and pairing gaps occur. We calculate these transitions 
by solving the Bogoliubov-deGennes equations both nu- 
merically and analytically using a degenerate shell pair- 
ing approximation for weak interactions and semiclassi- 
cal methods for many particles. An understanding of 
such elementary excitations is important if one wants to 
predict the response of the gas to various single particle 
probes. Moreover, these elementary excitations are also 
the building blocks, if one wants to calculate collective 
properties of the superfluid gas. 

The paper is organized as follows: In sec. ^ we set up 
the basic formalism used in the rest of the paper. We then 
in sec. 



III| discuss the weak pairing regime, where Cooper 



pairs are formed only between particles within the same 
harmonic oscillator shell. The qualitatively different limit 
of strong pairing is discussed in sec. IV, where Cooper 



pairs are formed between particles residing in different 
shells. We finally summarize our results in sec. 



II. BASIC EQUATIONS 

We consider a dilute gas of N fermionic atoms of mass 
TO in two hyperfine states \a =t, i) trapped by a spher- 
ically symmetric harmonic oscillator (h.o.) potential at 
zero temperature. The numbers N„ of atoms trapped in 
each hyperfine state are assumed to be equal such that 
N = 2Na- as this is the optimum situation for Cooper 
pairing. We assume that the atoms interact through an 
attractive s-wave scattering length, a < 0. The Hamilto- 
nian is then given by 
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It should be noted that to reproduce the correct 2-body 
scattering properties one needs to use a pseudo-potential 
6{r)drr instead of a simple delta-function for the interac- 
tion term in Eq.(|l|). This is important when the pairing 
correlations are to be calculated We will return to 
this point in sec. |l|. In a non-interacting system the 
Fermi energy is 



Ep = [np + 3/2)fiw ~ (3iV)^/^ hw , 



(2) 



where np is the h.o. quantum number at the Fermi sur- 
face. For a large number of particles np — (3A^)^/'^. 

Interactions result in a mean field potential which is for 
a large number of particles given by the Thomas-Fermi 
approximation (TF) 



U{r) = 27r p{r) , 

TO 

where the particle density is 

p{r) = 4(r)/37r2 ~ po (l - ^V^^If - U{r)/Ep) 



3/2 



(3) 



(4) 



inside the cloud r 



< Rtf — aosc\/2np + 3, where 
^Jhjnujj is the oscillator length, and po — 
{2npY/'^ /Zir'^a^gsc is the central density @]. kpir) is the 
TF Fermi wave-number and the Fermi energy is 



Bp = h^kl{r)/2m + (I/2)mw^r^ + U{r) 



2„2 



(5) 



In dilute systems the interactions contribute a mean field 
weak compared to the confining potential, i.e., U{r) can 
be neglected in Eq. (^. 

The system is dilute when the Fermi energy is large 
compared to the mean field energy, Ep ^ [/, or equiva- 
lently kp\a\ ~ /9^/'^|a| ~ ny^\a\/aosc <C 1. Dense Fermi 



I 



systems, kp\a\^l, are studied in Refs. |||. In very di- 
lute systems also the h.o. energy exceeds the mean field 

potentials, huj ^ U, or equivalently n^J^\a\/aosc 1- 

The h.o. levels are highly degenerate with states hav- 
ing angular momenta I = n, n — 2 . . . , (n even) or 
I — n,n ~ 2 . . . ,1 {n odd) due to the U (3) symmetry of 
the 3D spherically symmetric h.o. potential. However, 
interactions split this degeneracy. The splitting can be 
calculated perturbatively in the very dilute limit. From 
the radial h.o. wave-function Ti-npiif) for the state with 
angular momentum I and {np — 0/2 radial nodes in the 
h.o. shell np, the single particle energies with respect to 
the Fermi energy are within the WKB approximation 0] 



Cn^.i = / U{r)\'Rni{r)\'r'dr 



4^/2 
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The mean field energies are proportional to the coupling 
and vt^p"^ and are split like rotational bands with a pref- 
actor that is smaller than the average mean field energy 
by a factor ~ 3/16. The relative small splitting reflects 
the fact that the mean field potential is almost quadratic 
in r and the anharmonic terms therefore small. 

We plot in figure 1 the different regions of pairing for 
the trapped g function of the number of particles 

and the coupling strength. 
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FIG. 1. Tlie different level structures and pairing gaps in 
traps are displayed vs. the two dimensionless parameters: 
The energy quantum number nF of the highest shell occu- 
pied and the interaction strength \a\/aosc- The dotted line 
separates the dense and dilute trapped gases. Above the blue 
full line A = fttj, the system has coherence length smaller than 
the system size and pairing occurs between particles in differ- 
ent harmonic oscillator shells. Below, a region exists where 
there is only pairing between particles in the same shell but 
many h.o. shells contribute to the total pairing field. Below 
the purple line Gln(nF) — Q.lhw, there is only pairing in the 
shell at the chemical potential leading to the gap A ~ G for 
small systems whereas pairing takes place in a single /-level 
only for large systems. It should be noted that the transitions 
between the different regions are smooth. 



We shall below calculate pairing energies by solving 
the Bogoliubov-de Gennes (BdG) equations numerically 
and also analytically in certain limits exploiting known 
properties of the known h.o. wave- functions and single 
particle levels. The BdG equations read ||ict| : 



£;^u^(r) = [Ho + U{v)]u^{y) + A(r)w^(r) 
E^v^{v) = -[Hq + U{v)]v^{r) + A(r)u^(r) 



(7) 



Here, Ho = -{h^ /2m)\l'^ + C/o(r) - /ip, U{r) = 
.g(^/'J(r)i^cr(r)) — gPa{^) is the Hartree potential, and 
i/jai^) is the field operator that annihilates an atom 
in spin state a at position r. The coupling constant 
is g = Anah^ /ra with a < for an attractive inter- 
action leading to pairing, and ^p denotes the chemi- 
cal potential. The pairing field is defined as A(R) = 
— ghmr^o (R -|- r/2)?/;|(R — r/2))] and can be 

taken to be real 0]. The quasi particles with energies 
are described by the Bogoliubov wave functions w,, (r) 
and w,j(r). 

For spherical symmetry, the Bogoliubov wave functions 
from Eq.(|^) have definite angular momentum quantum 
numbers Im and can be written in the form u,-iim{'^) = 
r^^Ur,i{r)Yimi0,(l)) and u^/„i(r) = r^^Vr,i{r)Yi„r {0, (f)) 
with an energy energy i?^/. To obtain analytical results 
for the pairing energy in some limits, it is convenient to 
expand the Bogoliubov wavefunctions in the complete set 
of unperturbed spherical Harmonic oscillator functions 
by writing u,,/„i(r) = I]„ M7)n/</'°im(r) and -y,,im(r) = 

Y.„^r,nl(p°nimi^) with 0„;m (r) = TZnl{r)Yim{0 , (j)) . Wc 

haveT^o'/'VW = Cn<i™(r) with 6, = {n + 3/2)haj- ^p 
[I = n,n — 2 . . . ,0 (n even) oi I — n,n — 2 .. . ,1 (n odd)] . 
Using this expansion, the BdG equations become for a 
given Im: 



EniUnni = 'Y^[{nl\H\n'l)ur,n'i + {nl\/:^{r)\n'l)vr,n'i\ 
E,,iv,jni = ^[{nl\A{r)\n'l)ujjn'i - {nl\H\n'l)vnn'i\- (8) 



Here Ti = Hq + U (r) and the matrix elements are defined 
as {nl\A\n'l) = nni{r)Ann'i{r)r^dr. 
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We will examine certain regimes for which one can 
derive approximate analytical results from Eq.(^ [or 
Eq.(^] and thus for the Cooper pairing in a harmonic 
trap. We keep the temperature T = for the remaining 
part of this paper in order to examine the ground state 
properties of the gas. The nature of the paired state de- 
pends on the relative magnitude of the Hartree field U{r), 
the trap level spacing hu), and the pairing field A(r). 



III. INTRA SHELL PAIRING 



for the pairing field becomes: 

n,l 

where the factor 21 + 1 comes from the summation over 
the angular quantum number to. The nature of the self- 
consistent solution of the BdG equations depends on the 
position of the chemical potential /i p relative to the trap 
levels. We will now examine two opposite cases of inter- 
est. 



In this section, we examine the case when the interac- 
tion is so weak, that it only induces pairing between par- 
ticles within the same harmonic oscillator shell [n — n' 
in Eq.(H)], hence the term intra shell pairing. In this 
regime the pairing energy is smaller than the trap level 
spacing 2hLU between states with the same angular mo- 
mentum quantum number I [The Cooper pairing is be- 
tween states with the angular quantum numbers {l,m) 
and (Z,— to)]. We ignore for the time being the split- 
ting of the harmonic oscillator shells due to the Hartree 
field. The harmonic oscillator shells therefore have the 
energies ~ ^„ = (n -I- 3/2)huj — fj,p with the angular 
momentum quantum number Z = n, n — 2 . . . , (n even) 
or Z = n, n — 2 . . . , 1 (n odd) . The strength of the pairing 
between states with quantum numbers {n,l) and {n',l) 
in Eq.(H) is given by the matrix element 



{nl\A{r)\Ti'l) 



drr^A{r)nniir)nn'i{r). (9) 



The radial functions 7?.„;(r) in general change sign as a 
function of r. In the ground state however, A(r) has a 
definite sign (chosen to be positive in the present case) 
as any oscillations in the pairing field in general costs 
energy. We therefore see that the "inter-shell" matrix 
elements with n ^ n! will be suppressed due to the oscil- 
lating phase of 7?.„;(r )??.„' ;(r) as compared to the "intra- 
shell" matrix elements with n = n' . The weak pairing 
regime considered in this section is characterized by the 
fact that the inter-shell matrix elements can be ignored. 
A simple perturbative calculation reveals that the param- 
eter determining the importance of the inter-shell matrix 
elements is \{nl\A{r)\n'l)\/fiuj. When this parameter is 
small these matrix elements can be ignored and we are 
in the intra-shell pairing regime. For stronger coupling, 
this parameter is no longer small and the inter-shell ma- 
trix elements changes the solution to the self-consistent 
and thus non-linear gap equations qualitatively. We will 
describe this in detail in sec. [V. 

Ignoring all inter-shell matrix elements, Eq.(^) splits 
into simple 2x2 matrix equations for each pair of quan- 
tum numbers (n, I). For a given set of quantum numbers 
(n, Z), we obtain the usual solution Uni = (1 + (,n/ Eni)/2, 
Vni = (1 - £.nlEni)/2 with E^i = (e^ + A2,)i/2 and 
Ani — {nl\A{r)\nl) . Ignoring for the moment the regu- 
larization of the gap function, the self-consistent equation 



A. Partly filled Shell 

When the chemical potential is at a harmonic oscil- 
lator shell, i.e. ^ip = {np + 3/2)huj, this shell is partly 
filled for T — 0. This is the typical situation physically, 
whereas the case of a completely filled shell correspond- 
ing to lip = {np + 2)hu! a.t T = occur only for "magic" 
numbers of trapped particles. In this section we will pro- 
vide a analytical solution to the pairing problem when 
the interaction is so weak that only intra shell pairing is 
significant and iip = (np + 3/2)hu;. We will show that 
for very weak interaction the main contribution to the 
pairing correlations originate from the partly filled shell 
at the chemical potential (n = np), whereas for stronger 
interactions the contribution to the pairing from shells 
away from the chemical potential {n ^ np) is dominant. 
Using Eq.(^) and Eq.(|o|), the self-consistent equation for 
the intra shell pairing matrix elements reads 



Anl = 

E 

n' F ,1' 



r drr^TZl,ir)n. 

It •'0 



2 



(r) 



21' + 1 



4tt 



in' 



drr^nUr)nl,,{r) 
(11) 



where the first term singles out the contribution from 
the shell at the chemical potential. Notice that for 
this shell, S^npi = such that Enp = A„^;. Equation 
(11) gives the strength of the pairing as a function of 
the angular momentum I of the particles forming the 
Cooper pair and their shell quantum number n. The 
Z-dependence of the pairing, determined through the in- 
tegral J drr^ A{r)TZ\i{r) , is rather weak; we will later 
show that \An,i=n — A„_i=o|/A„,i=o <C 1. It is therefore 
not important for the present analysis and we define the 
Z-averaged pairing strength as A„ = + l)A„i/r2„ 

with f]„ = Y.ii'^l + 1) = {n + l)(n + 2)/2 being the 
(single spin) degeneracy of the harmonic oscillator shell 
with energy (n -I- Z/2)hw. For n 1, we can use the 
Thomas-Fermi identity 



E 



(2^ + i; 

47r 



K'lir) 



dpair) V2n' + 3 



dn' 



l-ryR^^p{n') (12) 



with RTp{n') = y/n/ + i/2ao8c- Here, pa is the density 
from a single hyperfine state when shells up to energy 
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(n'+3/2)fia; are occupied; dpair)/dn' therefore should be 
understood as the derivative of the density with respect 
to the highest harmonic oscillator shell occupied. The 
pairing strength A„ depends weakly on n for the shells 
around n = np which mainly contribute to the pairing 
and we therefore approximate A„ ~ A„^ = A. Equation 
(11) then becomes 



n<2nF 



A 



A = G+i^G y ^= , 



where we have defined 

\g\ J drr^dnPair)\n^]^ _ 32V2nf + 3 



G 



2 / drr'^dnPa{r)\nf 



157r2 



-huj . 

= (14) 



The correction factor i' arises because the overlap inte- 
grals in the second term in Eq.(ll) depend on the shell 
number n' . We approximate this dependence by averag- 
ing the factor •\/n', which gives v = /q^"^ n^/'^dn/2ri?p'^ = 
2^/^/3. We have introduced the cut-off n < 2nF in the 
sums which models as a first approximation the more rig- 
orous regularization procedure described in Ref. that 
is required for a delta-function pseudo-potential. Such a 
procedure can for the present purpose be approximated 
by a cut-off oin < 2nF essentially due to the fact that the 
Lippman-Schwinger eqn. for two particle scattering ap- 
pearing in the regularization procedure contains a single 
particle Green's function which is odd around the Fermi 
surface whereas the contribution to the pairing is even 
due to the particle-hole symmetry. 

When A htu we can expand and find the gap 



G 



(15) 



Here, 7 ~ 0.577 is Eulers constant. For very weak in- 
teractions, Gln(nF) ^ ^1^, the gap is simply A = G. 
This corresponds to the limit where pairing only occurs 
in the shell right at the chemical potential (n = np)- 
The result A = G was also obtained in ref. M using a 
different approach (the seniority scheme). As G\n{np) 
increases the gap increases and A > G due to intra shell 
pairing in shells away from the chemical potential. For 
Gln(nF) ~ 1/2, Eq. ( |l5|) diverges. However, at the same 
time the shell pairing ansatz leading to Eq. (11) breaks 
down as will be described in detail in the following sec- 
tion. 
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Eq.(29) 
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Eq.(15) 
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n[exp(y)nji'^ 
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0.2 0.3 0.4 0.5 
G [In trap units] 

FIG. 2. The lowest quasi particle energy i5i=o divided 
by the effective coupling strength G as a function of 
G ln[exp(7)n_p] (in units of huj) (a). Figure (b) depicts E/huj 
as a function of G/hcu. hf = 51.5/ia;. Plusses and crosses are 
with and without the Hartree field respectively. 



In fig. (|^) , we plot the lowest I averaged quasi particle 
energy E ^ A (since — 0) divided by the effective 
coupling strength G as a function of Gln[exp(7)ni?]. The 
solid curve is calculated from Eq. ( p^ ) whereas the x 's are 
results obtained from numerical calculations as described 
in ref. We have chosen np — 50 for the numerical cal- 
culations corresponding to ^ 2 x 10^ particles trapped in 
each of the two hyperfine states. Note that according to 
Eq.(|l5|), E/G in terms of Gln[exp(7)ni;'] is a universal 
function independent of the size of the system (i.e. in- 
dependent of np) as long as the system is in the intra 
shell pairing regime. In this regime, we see that there 
is excellent agreement between Eq.(p^) and the numeri- 
cal results. For stronger coupling, the intra shell pairing 
assumption breaks down and the numerical results dif- 
fer qualitatively from Eq.(15). The nature of the Cooper 
pairs undergoes a transition to a region where the pair- 
ing between particles in different harmonic shells is sig- 
nificant. The transition in the nature of the pairing is 
clearly visible in fig. (^) since E/G changes from being 
an increasing to a decreasing function of Gln[exp(7)rii?] 
when the intershell pairing becomes dominant. We will 



examine the intershell pairing in detail in section IV. Of 
course, the pairing strength is a monotonic increasing 
function of the coupling strength G as shown in fig. (||b). 
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For weak pairing, [Gln(n_F) <C fko] we find A = G in 
agreement with the analysis above. In this regime there 
is only significant pairing in the shell at the chemical po- 
tential. For stronger interaction with Glii{np) ^ O.lhiu, 
the Cooper pairing is still only between particles within 
the same shell, but shells away from the chemical po- 
tential (n np) contribute significantly to the pairing 
energy as predicted by Eq.jTsI). 

The reason that shells away from the chemical poten- 
tial contribute to the pairing for stronger coupling is 
that it becomes energetically favorable to partially oc- 
cupy these shells in order to make a coherent state with 
stronger spatial overlap between the particles forming the 
Cooper pairs. 

We will now briefly examine the dependence of the 
pairing strength A„y; given by Eq.(ll) on the angular 
momentum quantum number I. It is assumed that we 
are in the very weak coupling regime [ln(nF)G ^ 1] 
such that such that only the shell at the chemical po- 
tential n — np contributes to the pairing. Using Eq.(12), 

T^nFi^npi''') ^ 6xp(~''^/2«Lc)! the semiclassical 
approximation rTZnpi=Q{r) ~ [k(r)]~^/'^ sin(J''^ dr'k{r') + 

(/>) with k{r) = \/2np — /a^sc^ the radial integral in 



Eq.(ll) can easily be solved for / 
I — np yielding 



A„ 



G 



G 



15/16 
157r/(32V5) 



{np 

I 
I 





np 



and 



(16) 



with G given by Eq.(0. From Eq.(|l^), we see that 
|A„y,;=„p - Anp,i=o\/A ~ 0.10 < 1 thereby justifying 
the / averaging of the pairing strength described above. 
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FIG. 3. The lowest quasi particle energies Enpi — An^i (in 
units of G) as a function of the angular momentum l/np for 
= 51.5huj and G = 1.75 x IQ-^huj. 

Figure (|^) shows the /-dependence of the lowest quasi 
particle excitation energy En pi = An pi for the same 
parameters as for fig. (||) and G/Hlj = 1.75 x 10"'' <C 
\n{np)~^. The x's indicate the numerical results, the *'s 
indicate the results given by Eq.(p^. We also plot the 
I averaged value A = G as a solid line. The numerical 
results agree very well with the analytical predictions of 
Eq. (^6|) up to 1/np corrections. We notice that the 
/-dependence is weak thereby justifying the /-averaging 
of the pairing assumed above. Again, one should note 



that fig. (g) is universal in the sense that it is indepen- 
dent of the size of the system given by np; one always 
has < l/np < 1 and Enpi/G is independent of np for 
\n{np)G <C 1 as can be seen from Eq. 



B. Filled Shell 



We now examine the case when the chemical poten- 
tial is exactly between two harmonic oscillator shells, i.e. 
pLp = (np + 2)huj. This case corresponds to shells up to 
and including n — np being completely filled and shells 
with n > np empty when there is no pairing. As we 
will see, this situat ion d iffers qualitatively from the case 
considered in sec. ( Ill A ) in the sense that the pairing is 
zero below a certain critical coupling strength Gc- Again, 
we are in this section in the weak pairing regime such 
that only intra shell pairing is significant. Thus, Eq.(p^) 
is still valid and leads to a gap equation analogous to 
Eq.(ll), however, excluding the first term corresponding 
to pairing in the shell at the chemical potential. It is 
the absence of this term that leads to the vanishing of 
pairing below a certain finite critical coupling strength. 
Following steps completely analogous to the ones leading 
from Eq.(ll) to Eq.(|l|), we obtain the gap equation 



, 2n 

nui X ^ 

IJg^ ^ 



1 



np 



1/2)2 + (A/fi^)2 ^^^^ 



with G still defined by Eq.(^4|) and the scaling factor 
1/ = 2^/2 y'3 again approximating the n dependence of the 
radial overlap integrals. The critical coupling strength Gc 
below which there is no pairing is found by setting A = 
in Eq.(p^). Neglecting terms of 0{np^), we obtain 



G,= 



25/2 \-a{4:e^np) ' 



(18) 



Using Eq.(p^ and expanding Eq.(p7|) in the small pa- 
rameter A/huj, we find 



A 



G <Gr 



23/275(3) 



1/2 



with C(3) = EZi 



1.202. Equation (h9h should 



be compared to Eq.(15). We see that the main difference 
is that when there is a harmonic oscillator shell at the 
chemical potential, the pairing occurs as soon as there is 
an infinitesimal attraction between the particles. When 
the chemical potential is between two shells, there is only 
pairing when the coupling is above a certain critical cou- 
pling strength Gc as there is zero density of states at the 
chemical potential. 
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C. Effect of the Hartree field 




G/G 

c 

FIG. 4. The lowest I averaged quasi particle energy E/hu) 



(a) and the pairing strength A/ftoj 
as a function of G/Gc- 



In analogy with fig.(||), we plot in fig.(§) the lowest 
I averaged qua si p article energy as a function of G/Gc- 
The line is Eq.(|l9|) and the x 's are obtained from numer- 
ical calculations with fip — 52huj. As before, we conclude 
from Eq.(|l9|) that A as a function of G/Gc is indepen- 
dent of rii? in the intra-shell pairing regime. Of course, 
the size of the system is important when determining Gc 
as is apparent from Eq.(p^). Note that the lowest quasi 
particle energy E is calculated as = -I- A^)^/^ 

with S^np = fiuj/2 and A given by Eq. (|l9|) . To en- 
able a more direct comparison between the numerical re- 
sults for the pairing strength and Eq.(|l9|), we also plot 
A = {E'^-S.l )^/^ The observable in an experiment that 
probes single particle excitations is of course the excita- 
tion energy E and not A. We see that there is excellent 
agreement between the numerical results and Eq.([l8|)- 
( p^ ) for A/hiu < 1/2. The approximate analytical for- 
mulae predict very well the critical coupling strength and 
the size of the pairing for small A/h uj whe re only intra 
shell pairing is significant. As in sec.(IIIA), the pairing 
between particles in different harmonic oscillator shells 
become important when A/huj ^ 0{1) and the intra shell 
pairing theory breaks down. 

Pairing between just two shells was studied in the con- 
text of nuclear structure where similar results were found. 
Comparing to experimental data pd| only qualitative 
agreement was found. The nuclear spectra are compli- 
cated by stronger anharmonic mean field potentials and 
the finite range of the nucleon-nucleon interaction as well 
as spin and spin-orbit dependences. For atomic traps 
the results presented here should be more accurate as 
the pseudo potential approximation for the atom-atom 
interaction is expected to work very well for low atomic 
densities, and since the mean field is weaker. 



We have sections [II A- [II B ignored the Hartree field 

of 
In 



which induces an energy shift as given by Eq.(| 
the levels within a given harmonic oscillator shell 
this section, we will examine the effect of this mean field 
induced spread in the harmonic oscillator shell energies 
on the pairing. The gap equation including the Hartree 
field becomes in complete analogy with Eq.(ll) 



n'l' 



21' + 1 A„ 
47r 'Wr, 



drr'Rii{T)Ri,i,{r) 



with quasi particle energies 



E^^ 



nl 



rd 



(20) 



(21) 



and Rni{r) being the solution to the radial Schrodinger 
equation including the Hartree field. From Eq.(^, we see 
that the energy spread of a given shell due to the Hartree 
field is |^„i=o,i ~ ini=n\ rt'/'^a/ a^sc whereas the size of 
the pairing assuming completely flat shells is given by 
Eq.(|l5|) [or Eq.(|l9|)]. Thus, for large np the Hartree field 
splits the different / levels within the same shell suffi- 
ciently so that the pairing occurs only within one (or a 
few) energy levels (^npi close to the chemical potential 
Assuming that one level S^npi is at the chemical poten- 
tial (i.e. £,npi = 0), only this level will contribute to the 
pairing for weak coupling and Eq.(^) becomes 



A. 



21 + 1 
in 



(22) 



This equation is straightforward to solve as we simply 
have to evaluate the radial integral. For instance, for 
I = Hp, we obtain: 



Anpl—rif 

huj 



\a\/ao 



(23) 



Equation (23) agrees with the results derived in ref. |9[| 
using the seniority scheme. 

For Eq.(p3|) to be vahd, we need Anpi <C linpi ~ 
6if'±2|- Otherwise, the levels fnj^;±2±4... away from 
the chemical potential will contribute significantly to the 
pairing and the assumption that only particles in the 
level S^npi = form Cooper pairs leading to Eq.(23), 
will be incorrect. In the discussion following Eq.([5), 
we concluded that for ln(ni;')G > O.lhuj levels away 
from the chemical potential with normal phase energies 
— {n — nF)hiiJ contributed significantly to the pair- 
ing (see Fig. |2|). Equation (^ with one level at the 
chemical potential {£,npi = 0) is completely analogous to 
Eq.(ll); ignoring the contribution to the pairing from 
n ^ np shells in Eq.(pO|) the essential difference between 
the two equations is the level spacing between the normal 
phase energies which is huj in Eq.(ll) and ^npi — CnFi±2 
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in Eq.(pO[). Thus, a similar analysis on Eq.(gO|) yields 
that the contribution from states S^npi' ^ to the pair- 
ing can be ignored ifj^n^; - 6if;±2| > '^npi In(nF) with 
^npi given by Eq.(p3). For I = np which gives the max- 
imum l^npi — ^npi±2|, this condition yields [using Eq.(|^) 
and Eq.(p3|)] np ^ 200. That is, only for rather large 
systems with ^p ^ 200huj (corresponding to iV 10^ 
particles trapped), is the pairing for weak coupling only 
significant in the level ^nj,ij, = at the chemical poten- 
tial and Eq.(p3|) is valid. If the chemical potential is at 
a level with I < np (i.e. £,npi = with I < np), the 
Hartree splitting of the levels is smaller than for I = np 
and the condition for single {np,l) pairing becomes even 
more strict M. 



The assumption in sec. HI A -[II B that the Hartree field 



can be neglected such that the harmonic oscillator shells 
are completely flat, is consistent if A > |Ctif'=o,i ~ 
in^i=n^\ with A given by ^q.{^ [or Eq.(|l|)]. From 
Eq.(ll), we have |6yi=o,i - infi^nfl ^ (5/32)n_FG. Us- 
ing this and Eq.(|lq) in the very weak coupling regime 
with A = G, we get the condition np < 6 for the valid- 
ity of neglecting the Hartree field. For stronger coupling 
where pairing in the shells away from the chemical poten- 
tial is important, we of course obtain that the neglect of 
the Hartree field is consistent for larger np, since Eq.(|l5|) 
gives A > G. In particular, for A ~ 0.5huj which is the 
approximate maximum value for which the intra shell 
pairing assumption holds, we obtain that Eq.(p^) is valid 
for systems with np < 30. 

In order to examine the influence of the Hartree field, 
we have plotted in fig.(||)(b) as -l-'s the lowest quasi 
particle energy (for which = such that En^-i = 

\J~^i~i^r'^i = A„y/) including the effect of the Hartree 
field in order to compare to the results neglecting the 
Hartree field. We keep the number of particles trapped 
the same with or without the Hartree field. Since np — 
50 > 30, according to the analysis above the inclusion 
of the Hartree field should yield a significant modifica- 
tion o f the p airing as compared to the analysis presented 
As expected, we see that in the intra shell 



[HA 



pairing regime the Hartree field suppresses the pairing 
as it introduces a splitting of the levels in the shell at 
the chemical potential such that the density of states 
for ^jij = is reduced. For stronger interaction, where 
there is pairing between particles in different harmonic 
oscillator shells, the effect of the Hartree field is to in- 
crease the pairing strength. This is because in this regime 
the Hartree field compresses the gas increasing the den- 
sity p{r = 0) thereby increasing the effective interaction 
strength ^^(r = 0)|a| [See Eq.(|)]. We see that Eq.(|l|) 
gives a qualitatively correct description of the pairing 
strength for this set of parameters since rtj? = 50 is an 
intermediate value where the pairing in the intra shell 
regime is not in the limit where the Hartree field can be 
ignored, but also far away from the limit where there is 
only pairing within the level = 0. 



IV. PAIRING IN THE A(i? = 0) > ftoj LIMIT 

In this section, we will derive an analytical expression 
for the lowest single particle energy in the regime, where 
the pairing field A(r) is much larger than huj, yet still in 
the dilute limit so that A(r);Si?i?. We will show that the 
strong pairing between particles in different shells leads 
to a b reak-down of intra shell pairing ansatz of secs.III A- 
III B; we are in the inter shell pairing regime. 

When A(r = 0) ^ /iw the local or TF coherence length 
in the center of the trap 



^{r = 0) = 



kp{r = 0) 
7rmA(r = 0) 



(24) 



is smaller than the size of the system. The coherence 
length in general determines the length scale over which 
the pairing field can vary. In this regime, the shell struc- 
ture of the trap levels can be ignored and the gas behaves 
in many ways as a quasi-homogeneous bulk system. In 
nuclei gaps are always smaller than hw by an order of 
magnitude [|Tl[| . Thus trapped atoms provide a bridge be- 
tween condensed matter and nuclear physics for studying 
superfiuid properties. 

The pairing field A(r) in Eq.(|7|) in general increases 
with r ^ in the bulk regime. When A(r = 0) ^ Two, 
we have ^Bcsif = 0)/Rtf ^ 1- One can in this regime 
use a semiclassical theory developed in Ref. to ob- 
tain an analytical expression for the lowest quasi parti- 
cle energy. As discussed in Ref. |T^], the lowest exci- 
tations are concentrated in the region around the mini- 
mum of A(r) + Uoir) — fip. The lowest excitations have 
energies E^ii significantly lower than A(r = 0) and can 
be denoted "in gap excitations" . Their wave-functions 
are suppressed in the regions where E^ii < A(r) and 
Erji + fip < Uo{r). To illustrate this, we plot in Fig. 
(^) the four Bogoliubov wave functions [u,,;(r), i;^/(r)] 
with lowest energy E^ji for I — 0, ^p — 51.5/ia; and 
G ~ O.Mhuj. 




FIG. 5. The lowest 4 Bogoliubov wavefunctions Uni{r) and 
Vrjiir) in units of a~^^ for Z = 0, A(r), and Uo{r) — fip in units 
of hu!. 
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The solid lines are the Ujji=o{r) and the dashed lines 
are the Vrii=o{r). The wave functions are centered hori- 
zontally along a y- value equalling their excitation energy 
in units of huj. We also plot A(r) and Uo{r) — fip 
illustrate how the lowest quasi particle excitations with 
energies Ejj ^ A(r = 0) are centered at the edge of the 
cloud where A(r) + Uo{r) — fip has a minimum. 

For Ejj A{r = 0), a semiclassical analysis yields the 
following quantization condition for the I = states [n2[ : 



(77 



— )—fkj 
2' 2 



Rtf 



dr 



A(r)2 



\/ Rtf ^ 



(25) 



with 77 = 0, 1, ... . Here, Vc is the left turning point de- 
termined by A(rc) = E,j. For A(r = 0) ^ hoj we have 
l—Tc/ Rtf ^ 1- The centrifugal barrier, present for I > 
plays a minor role since the Bogoliubov wavefunctions are 
suppressed by A(r) in the region where the centrifugal 
potential is significant. Therefore the /-dependence of E^i 
is weak and only when I becomes comparable to tif does 
the centrifugal potential affect the lowest energy states 
for a given I by increasing its energy. 
In a uniform system the pairing gap is 



A 



2m 



(26) 



The prefactor k = 8e~^ within the usual BCS theory but 
is further reduced by a factor (4e)^^/'^, when induced 
interactions are included jl^ . 

In a non-uniform system such as the h.o. traps we 
can use the local density (or TF) approximation in the 
A(r = 0) ^ huj regime, where the discrete level structure 
is of minor importance as compared to the strong pairing 
field. A Thomas-Fermi treatment yields 



A(r) ~ K^lF{l-r^/RTF)e 



-7r/2fcF(r)|a| 



(27) 



The gap decreases rapidly near the surface and we can 
therefore neglect the gap in the region Tc < r < Rt f [see 
Fig. (||)]. The semiclassical quantization condition thus 
simply gives 



vr 

E„ — — 



2 a,TCCos{rc/ Rtf) 



(28) 



Here, the turning point given from Erj = A(rc) is now 
determined from the gap of Eq. (^7|) . We thus obtain in 
the A{r = Q) ^ hu regime 



G 



•i2E„ 



157r2(7?-h 1/2) 



In 



k(?7 



l/2)27r2 HFihiof]]^ 



E^ 

V 



(29) 



from which the quasi particle particle energy can be 
found by inversion from the coupling strength G defined 
in Eq.(g). 

The lowest quasi particle energy, Eq, from Eq. ( p9| ) 
(plotted as a dashed line) is compared with the numerical 



results in fig.(0) and fig.(|j). We see that the agreement 
between the semi classical prediction and the numerical 
results is good for E ^ huj where A(r = 0) ^ hw. Note 
that in the inter shell pairing regime, the position of the 
chemical potential with respect to the Harmonic oscil- 
lator levels is unimportant as the lowest quasi particle 
energy is determined by the geometry of the minimum of 
A(r) -I- Uo{r). Thus, Eq.(p9|) agrees well with the numer- 
ical results for the parameters relevant for both fig.(|^) 
(partly filled shell) and fig.(§) (filled sheU). Combined 
with the theory presented in sec. (HI), we therefore have 
a theory describing the strength of the pairing fairly well 
in the regimes A < ftu;/2 and A{r — 0) ^ huj. In the 
weak coupling regime, the pairing is between particles 
within the same shell whereas for stronger interaction, 
there is pairing between particles in different shells and 
the lowest energy excitations are states localized in the 
surface region around the minimum of A(r) -I- Uo(r). 

Equation (29) for the lowest quasi particle energy 
should be compared with Eq.([27|) giving the size of the 
pairing field in the center of the trap. We see that con- 
trary to the homogeneous case, the lowest pairing energy 
in general is different from the size of the pairing field 
even in the quasi homogeneous bulk limit; it grows much 
slower with G than A(r = 0). 



V. CONCLUSION 

In this paper, we considered the properties of the ele- 
mentary excitations of a superfiuid trapped atomic Fermi 
gas. The nature of the paired state was shown to un- 
dergo transitions between several different regimes as 
the number of particles (or the h.o. quantum number 
uf — {3Ny^^) and the strength of the interaction are 
varied. Three regions emerge: 

• For weak coupling, G ln(rti?)^0.1, the Cooper pairs 
were shown to be formed between states in the h.o. 
shell at the Fermi surface only. In this limit, it was 
demonstrated that both the position of the chemi- 
cal potential with respect to the harmonic oscillator 
shells and the size of the Hartree field, are strongly 
infiuencing the nature of the paired stated. The 
pairing gap is A = G for sufficiently few particles 
in the trap, uf^W, but for more particles the mean 
field splits the single particle levels and reduce the 
pairing. For Uf ^ 100 the pairing is reduced to 
that in a single /-level. 

• For interactions such that 0.1^Gln(ni?)/?iijj^l/2 
Cooper pairs are still essentially only formed be- 
tween states within the same shell. However, there 
is pairing in many shells besides that at the Fermi 
level yielding a gap significantly larger than A = G. 

• For strong interactions, Gln(nF)/Sw^l/2, the 
pairing field exceeds A(r = 0)}ijhuj. The coher- 
ence length is shorter than the system size, and 
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the gas can for many purposes be regarded as a 
quasi-homogeneous bulk system The Cooper pairs 
are formed between particles in different shells and 
the elementary excitations are confined spatially to 
the surface of the gas. 

It should be noted that the transitions between the re- 
gions described above are smooth contrary to thermody- 
namic phase transitions. The analytical formulae appro- 
priate for the various regimes were compared with nu- 
merical calculations showing excellent agreement. The 
results presented in this paper suggest, that by varying 
the number of particles and/or their interaction strength 
the trapped atomic gases can behave either essentially as 



giant nuclei (Sec.^) or as bulk systems (Sec.lV). They 
thus provide an intriguing link between uniform systems 
explored in condensed matter physics and finite systems 
relevant for nuclear physics, quantum dots, etc. 
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